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Abstract 

We consider the linear growth of matter perturbations on low redshifts 
in some f{R) dark energy (DE) models. We discuss the definition of dark 
energy (DE) in these models and show the differences with scalar-tensor DE 
models. For the f{R) model recently proposed by Starobinsky we show that 
the growth parameter 70 = 7(2: = 0) takes the value 70 ^ 0.4 for Om.o = 0.32 
and 7o — 0.43 for r^m.o = 0.23, allowing for a clear distinction from ACDM. 
Though a scale-dependence appears in the growth of perturbations on higher 
redshifts, we find no dispersion for 7(2;) on low redshifts up to z ~ 0.3, 7(2) 
is also quasi-linear in this interval. At redshift z = 0.5, the dispersion is still 
small with A7 ~ 0.01. As for some scalar-tensor models, we find here too a 
large value for 7^ = ^{z = 0), 7^ ~ -0.25 for Qmfl = 0.32 and 7^ ~ -0.18 for 
^m,o = 0.23. These values are largely outside the range found for DE models 
in General Relativity (GR). This clear signature provides a powerful constraint 
on these models. 
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1 Introduction 



The present stage of accelerated expansion of the universe [I] is a major challenge for 
cosmology. There are basically two main roads one can take: either to introduce a 
new (approximately) smooth component or to modify the laws of gravity on cosmic 
scales. In the first class of models, the new component dubbed dark energy (DE) 
must have a sufficiently negative pressure in order to induce a late-time accelerated 
expansion. One typically adds a new isotropic perfect fluid with negative pressure 
able to accelerate the expansion. In the second class of models one is trying to explain 
the accelerated expansion by modifying gravity, after all the universe expansion is a 
large-scale gravitational effect. It is far from clear at the present time which class of 
DE models will finally prevail and one must study carefully all possibilities. While 
all models are called DE models [2], the models of the second class are called more 
specifically modified gravity DE models. 

It is clear that in both classes we have many models able to reproduce a late-time 
accelerated expansion in agreement with the present data probing the background 
expansion like luminosity distances from SNIa. Probably this will remain the case 
even with more accurate data. However the evolution of matter perturbations will be 
affected differently depending on the class of models we are dealing with. This can be 
used for a consistency check in order to find out whether a DE model is inside General 
Relativity (GR) or not [3]. Therefore the evolution of matter perturbations seems 
to be a powerful tool to discriminate between models inside and outside GR In 
particular, the growth of matter perturbations on small redshifts seems a promising 
tool and it has been the subject of intensive investigations recently [5]. 

A particular class of modified gravity DE models are f{R) models where R is 
replaced by some function f{R) in the gravitational Lagrangian density [B] (for a 
recent review, see [7]). The idea of producing an accelerated stage in such models 
was first successfully implemented in Starobinsky's infiationary model [8]. In some 
sense it was rather natural to try explaining the late-time accelerated expansion by the 
same kind of mechanism, a modification of gravity but now at much lower energies. 
Some problems arising in these f{R) DE models were quickly pointed out, related 
to solar-system constraints [9] and instabilities ^lU]. It was then found that a very 
serious problem arises where it was not expected: many of these models are unable 
to produce the late-time acceleration together with a viable cosmic expansion history 
[llj . It was shown that many popular f{R) models, like those containing a R~^ term 
[12], are unable to account for a viable expansion history because of the absence of a 
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standard matter-dominated stage a oc is which is replaced instead by the behaviour 
a (xt2 [T3]. 

However some interesting cosmological models remain viable like the model re- 
cently suggested by Starobinsky [13], (for another interesting viable model see also 
|15j). It is interesting in the first place because of its non-trivial ability to allow 
for a standard matter-dominated stage. This goes together with the appearance of 
large oscillations in the past and the overproduction of massive scalar particles in the 
very early universe which, as noted already in [T3], can be a serious problem of all 
viable f{R) DE models. Another interesting property is the scale dependence of the 
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growth of matter perturbations deep inside the Hubble radius [IS],[IZ]- This scale 
dependence was actually used in [H] in order to constrain one of the parameters of 
the model. Of course this model also satisfies the local gravity constraints for certain 
parameter values. In the end a window in parameter space remains where the model 
is in principle physically acceptable. 

In this letter we want to assess the viability of this model with respect to the 
growth of matter perturbations (see e.g. [18] for other possible constraints and ap- 
proaches). Some viable f{R) DE models have their free parameters constrained in 
such a way that they cannot be distinguished observationally from ACDM. The situa- 
tion is different in this case. Our results show that the growth of matter perturbations 
on low redshifts z < 1 provides a discriminating signature of these models, able to 
clearly differentiate it from ACDM and also from any other DE model inside GR. 
These results confirm that the growth of matter perturbations can be used efficiently 
to track the nature of DE models. 



2 f{R) modified gravity cosmology 

We consider now f{R) models and discuss some general properties in connection with 
their observational viability. Many f{R) models have the surprising property that 
they cannot account simultaneously for a viable cosmic expansion history and a late- 
time accelerated expansion. This severely constraints viable f{R) modified gravity 
DE models. We consider a universe where gravity and the content of the universe are 
described by the following action 



S = / dt^x\/—g 



167rG, 



For the time being, is a bare gravitational constant, its connection with obser- 
vations depends on the theory under consideration. We concentrate on spatially fiat 
Friedman-Lemaitre- Robertson- Walker (FLRW) universes with a time-dependent scale 
factor a{t) and a metric 

ds^ = -dt" + a^{t) di^^ . (2) 
For this metric the Ricci scalar R is given by 

R = q(2H^ + H^ , (3) 

where H = ^ is the Hubble rate while a dot stands for a derivative with respect to 
the cosmic time t. DE models have a > on low redshifts. The following equations 
are obtained 

3FH^ = SnG, (p„ + prad) + ^{FR~f)~ 3HF , (4) 
-2FH = StvG, (p^ + ^prad^ + F-HF , (5) 
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where 

In standard Einstein gravity (/ = R) one has F = 1. The densities pm and prad satisfy 
the usual conservation equations pi = —3H{1 + Wi)pi with Wm = and Wrad = |- 

We summarize briefly the way in which the dynamics of f{R) models can be ana- 
lyzed. For a general f{R) model we introduce the following (dimensionless) variables 



X2 



f 



-3 = 6§5 = ^ + 2. (9) 



-4 = (10) 



From Eq. (jlj) we have the algebraic identity 



= ^3^]j-^"' = 1 - a;i - X2 - X3 - X4 . (11) 



It is then possible to write down the following autonomous system 

-1 - X3 - 3X2 + x\- XiXz + Xa , (12) 

X2(2x3 - 4 - xi) , (13) 

-2x3(x3-2), (14) 



dxi 
dN 

dX2 X1X3 

dN m 

dx^ xiX'i 



dN m 
dx4 



-- -2x3X4 + 3:1X4, (15) 
dN 

where = Ina. In eqs( !T3l[T^ . the quantity m corresponds to 

d\ogF RF' 

(a prime stands for derivative with respect to R) hence m is a priori a function of R. 
However it is easy to see that R can in turn be expressed in function of the variables 
of our autonomous system. Indeed we have 

^3 _ RF _ dlogf 

X2 / d\ogR ■ ^ ^ 

Inverting f|T7|) we can in principle express i? as a function of || = r and so m becomes 
in turn a function of r, m = m{r) which closes our system. 



3 



It was found in [13] that when radiation is neghgible (x4 = 0) there exists a critical 
point Xi = 0,X2 = — |,X3 = i corresponding to an exact matter phase 



^eff = -l-|^ = ^(l-2^3)=0 (18) 
= I. (19) 

This is the critical point called P5(m = 0) in [13j, it satisfies r = — 1. 

Let us comment about the physical meaning of these conditions. For r = — 1 or 
X2 + = 0, eq.(^ reduces to 

3FH^ = SttG, pm - SHF . (20) 

2 

Then we clearly have an exact matter phase (a ~ ta) for F = or xi = 0. We note 
further from eqs. fll3lll4p that X2 + 2:3 remains zero for Xi = 0. 

The condition m = is more subtle. From eq. (HM a stable matter phase requires 

Xi = 3 m . (21) 

It is straightforward to generalize (pT!) for arbitrary scaling behaviour of the back- 
ground 

xi = 3 m(l + Wes) ■ (22) 

Actually neither xi nor m can be exactly zero during a matter phase (which corre- 
sponds to a critical point of the system). This would require F' = 0, which reduces 
to General Relativity plus a cosmological constant. The system can only come in 
the vicinity of P5 and trajectories with an acceptable matter phase have Xi ~ and 
m ~ with xi ~ 3 m. In addition, from F' > both xi and m should be positive 
during the matter phase. Finally, it is easy to check that fl2T|) implies 

Roca-^ . (23) 

The situation is very different in scalar-tensor models where a standard matter phase 
is possible with F ^ because the energy density associated with the dilaton can 
mimic the behaviour of dustlike matter for a particular type of potential U{(j)) jl9| . 
It is this property which is used in [20]. More explicitly, for the Lagrangian density 

C = (f(<I>) R - g^''d^<i>d,<^> - 2/7(<l>)) + L^g^,) , (24) 

we have the Friedmann equations (in fiat space) 

$2 

3FH^ = SttG, + y + - , (25) 
-2FH = StiG, pm + ^'^ + F - HF . (26) 

It is possible to have a solution with F = and ^ + U oc by a suitable choice 
of the potential U [19]. Because of this property it is possible in the matter phase 
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to have a nonvanishing DE energy density Pde- The definition oi pnE requires some 
care. If we want pde^o obey the usual energy conservation for perfect fiuids, certainly 
a desirable property, then the choice is not unique. The Friedmann equations (H1I5!) 
and also fl25f26p can be recast in the form 

3AH^ = SttG, ipm + PDE{A)) (27) 
-2AH = 8nG. {pm + Pde{A) + pde{A)) . (28) 

where A is some arbitrary constant. We have written pde{A) (and pde{A)) as these 
quantities depend on the choice of A. For example from fHfSj) . neglecting radiation, 
we obtain 

8nG,pDEiA) = ^(FR- f)-3HF + 3iA- F) H'' (29) 
SttG, (pDEiA) + pDEiA)) = F-HF-2{A-F) H (30) 

Obviously, whatever choice is adopted the corresponding DE component obeys the 
usual conservation law of an isotropic perfect fiuid. For the representation (1271) . it is 
natural to introduce the cosmic relative densities 

3AH^ • ^^^^ 

A natural choice, especially if we are willing to compare modified gravity theories 
with General Relativity (GR), is to take 

^^G^, (32) 

to very high accuracy, where Gat is Newton's constant found in textbooks. In GR, it is 
the bare gravitational constant appearing in the action and it corresponds to the value 
of the gravitational constant obtained in a Cavendish-type experiment measuring the 
Newtonian force between two close test masses. As is well-known, this identification 
no longer holds outside GR. For modified gravity theories it is the quantity Gefr, a 
priori different from G*, which plays the role of gravitational constant and which 
is measured in a Cavendish-type experiment and clearly all acceptable models must 
have the same Ges today to very high precision. 

For scalar-tensor models we can make the choice A = Fq [19] as ^ ~ Gefr.o to 
very high accuracy because ujbd,o > 4 x 10^ from solar system constraints. As we will 
see later, in our f{R) model a Cavendish-type experiment corresponds to the limit 
Geflf = G^ (see [HI and eq. (1501) below). It is therefore natural to take ^4 = 1 in such 
a model. 

As noted in [T^ a consequence of the arbitrariness in the choice of A is that it 
introduces an arbitrariness in the sign of pde as well. In a model where Fq > F in 
the past, Pde remains positive if we take A = Fq. In our f{R) models however, we 
have Fq < F{z > 0) in an expanding universe with F — 1 at high redshifts. Hence 
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Pde{A) defined from fl271) can become negative in tlie past if we take A = Fq. This 
problem is avoided if we define pde from fl27|28l) putting A = 1 



3H^ 



SttG, {pm + Pde + Pde) 



(33) 
(34) 



~2H 



Tlie corresponding relative densities become VLi = \^2^ and actually ri^i^; — > in 
the past. Note that though it is always possible to write the Friedmann equations in 
the form fl33|34l) this should not conceal the fact that we are dealing with a genuinely 
modified gravity theory obeying eqs. (Hll5l) . The difference is absorbed in the micro- 
scopic definitions of the DE component. This shows again the need to go beyond the 
background expansion in order to discriminate modified gravity DE models from GR. 

3 Starobinsky's model 

It was shown that a viable cosmic expansion is rather problematic in f{R) models. 
As already mentioned in the Introduction problems related to instabilities [TU] and 
solar-system constraints [9j were pointed out soon after the first attempts to produce 
f{B) DE models. An even more unexpected problem was found later related to the 
background expansion as many f{R) models, including popular ones containing an 
inverse power of i?, are unable to reproduce a viable expansion history containing 
a standard matter-dominated stage [H]. The viability of these models was then 
systematically studied in [13]. All this shows that the construction of a viable f{R) 
DE model is something highly non-trivial. Still there are some f{R) DE models left 
that can account for a viable expansion history and still depart from ACDM. One 
such interesting model was suggested recently by Starobinsky 



with A, n > 0. We see that f{R) — > R for — > in contrast to models containing an 
inverse power of R, in addition in flat space the correction term to Einstein gravity 
disappears. The quantity Rc has the dimension of R and is a free parameter of the 
model, it corresponds essentially to the present cosmic value of R. Recently problems 
related to the high-curvature regime of this model were pointed out [2T] but it seems 
this could be cured by the addition of a term ^ as noted in IT5I. However the mass 
m is required to be very high with m ^ 1 GeV hence this term is negligible in the 
low curvature regime. In the present work we will show that it is the low curvature 
limit of this model, through the growth of matter perturbations on small redshifts, 
that could severely constrain this model. It is convenient to introduce the reduced 
dimensionless curvature R=^. We have for this model 




(35) 




(36) 
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This is hardly inversible but fortunately the situation improves for R ^ 1. In this 
regime we have 



f{R) 



R-XRc = R- 2A(cx)) 
1 



(37) 
(38) 
(39) 



F 



m{r) 



r 



-1 - XR-^ 
2n{2n+ 1) 




2n(2n+ 1)^ 



A 



0, m > (40) 



3F' 3mF 



(41) 



It is clear that when f{R) satisfies fl37|) . X2 + comes close to zero as FR — f — 
2A(cxo) <C 6H^. We see that M^, the scalaron mass introduced in [S], becomes very 
large in the past. As noted in [13], one should find a mechanism to avoid becoming 
too large in the early universe. We get also strong oscillations 6R in the past and we 
illustrate this with the oscillations induced in the behaviour of Wes on high redshifts 
already during the matter phase as displayed on Figure la. These oscillations could 
lead to overproduction of scalarons in the very early universe [H]. While the cosmic 
background today corresponds to i? ~ 1, the regime R^ 1, eqs. (P7l) - (HTi) . is quickly 
reached as we go back in time. In this regime r ^ — 1 and m ^ so that a viable 
matter phase is possible followed by a late-time accelerated stage. We have also in 
this regime F ^ 1 to very high accuracy. One finds that in this model accelerated 
expansion starts at Za — (0.5 — 1) depending on the value of flm,o, in agreement with 
observations [22j . 

In order to be viable, this theory should not be in conflict with gravity constraints. 
On scales satisfying today k <^ 2ttM, one will not feel the presence of a "fifth force" 
caused by the scalar part of the gravitational interaction. Hence provided the scalaron 
mass M corresponding to some experiment is large enough, the "fifth force" will only 
be felt on scales that are much smaller than those probed by this experiment. In 
particular one will not see it in a Cavendish-type experiment on scales larger than 
5 X 10~^cm provided n > 1, and this range can be increased to even lower scales 
for larger n [T3]. We refer the interested reader to the reference [H] for a thorough 
discussion of the properties of this model and of its potential problems. 

The evolution and growth of matter perturbations can also be seen as an exper- 
iment probing gravity. In this case the scalaron mass can be much smaller so that 
deviations from GR and the appearance of a "fifth force" can be felt on larger scales 
and even on cosmic scales. We will consider this in more detail in the next section. 

4 Linear growth of matter perturbations 

In a way analogous to scalar-tensor DE models, the equation governing the growth 
of matter perturbations on subhorizon scales is of the form [23j 



Pm 



m 



. 



(42) 
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Figure 1: a) On the left panel the effective (total) equation of state parameter Weff 
is shown for n = 2. We see that large oscillations appear in the past, already in 
the matter-dominated stage at rather low redshifts. b) On the right panel, the 
behaviour of the quantity ^"^^ is displayed. It is seen that some subhorizon scales 

satisfying k <^ aH in the past will gradually shift to the regime k ^ aH as the 
universe expands. In this asymptotic regime, the scalaron corresponding to the cosmic 
backgroung curvature is nearly massless and a fifth-force is felt on these cosmic scales. 
Note that today Hq 7 Hq. 



The expression for G^s in f{R) models is given below, see eq.(42). Many DE models 
outside GR have a modified equation of this type describing the growth of matter 
perturbations. For the models where it is valid, this equation includes all perturba- 
tions, the assumption is that we are in a regime deep inside the Hubble radius where 
the leading terms result in (H2l) . One can also have in principle more elaborate DE 
models yielding a different equation (see e.g. [21]). It is convenient to rewrite the 
corresponding equation satisfied by the quantity / = j|M^[25] 



where Qm is defined from (!3T|) with A = 1, and [16 



k^F' 



1 + 3^ 



The r.h.s. of is given by the expression AnGeff^ which obviously does not 
depend on (the constant) A so that the evolution of / does not depend on A either. 
If we write this r.h.s. using the definition (29) Vtm = ^s^j/a" with A 7^ 1, we will have 



2„l 



A-KGeff^ = \ A %f^f^m and we have in particular | A %f^m = | ||^^ ^m, with 



Qm defined in eq.(ll). When gravity is described by GR eq. fH3l) reduces to eq.(B7) 
given in [26] . 



^The reader should not confuse this quantity with the quantity /(-R), we beUeve such confusion 
will be easily avoided from the context where they appear. 
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A crucial point in this model is that G^s is scale dependent 



G,s = G,s{z,k) . (45) 



In other words, the driving force in eg. (14311 introduces a scale dependence in the 
growth of perturbations. The spatial variation on cosmic scales of Ge// increases 
the growth of matter perturbations as is seen very clearly from Figure 2a. As this 
increase is rather tightly constrained by observations, so will the model parameter n 
on which this increase depends [Hj. We will return in more details to this important 
point below. 

We can rewrite Ges in the suggestive way 

a 1 + 4- " 



Ges = ^ ^ , /If (46) 

(47) 




where 



H^ = — = 3FM^ = — . (48) 
F' m 

In equation (H3l) . G^s depends on the cosmic curvature R taken from ([3]). In the limit 
i? ^ 1 satisfied by the cosmic curvature R at high redshifts, we have 

1 



Geff = a ( 1 + ^ 1 ^» 1 • (49) 



Taking the (inverse) Fourier transform of fH9|) it is straightforward to recover the 
corresponding gravitational potential per unit mass V{r). Remembering that a po- 
tential oc - in real space yields a term in Fourier space, we recognize in f H9l) the 
gravitational potential in real space (per unit mass) 

G* / , , 1 -Mr 



V{r) = -^ ^l + -e-''^j . (50) 

The quantity H (and M{R)) can become small enough with the universe expansion, 
see Figure lb, so that some cosmic subhorizon scales can feel a significant fifth-force. 
As the universe expands, H is rapidly decreasing so that these deviations are felt in 
Poisson's equation on ever increasing scales. While in the past only scales very much 
smaller than those corresponding to cosmic scales today could feel deviations from 
GR, today Hq ~ Hq and deviations can be felt on essentially all subhorizon scales. 
We can distinguish two basic asymptotic regimes 

Geff = k:^aH, (51) 

3 F 

= % k<^aH . (52) 
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The quantity aH is rapidly increasing in the past. Hence some cosmic subhorizon 
scales {k ^ aH) will satisfy k -C aH in earlier times (large redshifts) and switch into 
the regime k ^ aH on lower redshifts as the universe expands, see Figure lb. For 
cosmic scales that are in this regime (15T!) . the scalaron mass appears negligible and a 
fifth force does appear which results in the factor |. 

An accurate evolution of / or 5m for arbitrary r.h.s. of eg. (1431) requires numerical 
calculations. But it is possible to find analytically the solutions in the two asymptotic 
regimes fl51|52l) during the matter stage when VLm ~ 1 and F 1. This regime is 
still valid until low redshifts. When C = ^^fim = constant we get constant growing 
mode solutions f = p with (see e.g. [20]) 



P 



^ (-1 + Vl + 24C) . (53) 



A similar result can also be obtained in the framework of chameleon models (see e.g. 
[27j ) and inside GR (see e.g. [28]). We have C ^ 1 for -C and C ^ | for 
k ^ aH, so that we obtain [14] 



5m a -C aH (54) 

5m OC a 4 k ^ aH . (55) 



For subhorizon scales that go from the first into the second regime, this can yield a 
scale-dependent increase (see Figure 2a) in the growth of matter perturbations. Note 
that this increase takes place on small cosmic scales. This increase in the growth of 
matter perturbations will induce a change of shape of the matter power spectrum 
P{k) inferred from galaxy surveys with a subsequent change of its spectral index nf^', 
P{k) oc k"-" . This can then be compared with the spectral index n^^^ derived from 
the CMB anisotropy data resulting in a possible discrepancy between both spectral 
indices. No significant discrepancy between the two values is allowed by present 
observations and we have the conservative bound (see e.g. |29j ) 

nf - < 0.05 (56) 

So this difference is rather tightly constrained and so is the model parameter n on 
which it depends. This results in the constraint n >2. Because increasing the model 
parameter n will decrease this discrepancy according to the analytical estimate [14] 



^gal_^CMB^ V33 5 

2(3n + 2) ' ^ ' 

which is accurate for a wide range of n values (only for large values one should resort 
to a more accurate numerical estimate), (I56f57p result in the constraint n > 2. 

Let us emphasize again at this point the differences with scalar-tensor DE models. 
In scalar-tensor DE models Ges does not depend on k (or on r in real space), at a 
given time it is the same Ges that enters the equation for the growth of perturbations 
and the gravitational constant measured in a Cavendish-type experiment. These 
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DE models have further a neghgible dilaton mass but they can comply with solar- 
system constraints if they satisfy ujbd,o > 4 x 10"^, where u!bd,o is the value today of 
the Brans-Dicke parameter, which yields ^ 1 [23], [30]. Indeed for scalar-tensor 
models one has for the Post-Newtonian parameter today Jpn = !^g^'o+2 solar- 
system constraints point at jpn ~ 1 close to its GR value jpn = 1- As f{R) models 
correspond to scalar-tensor models with ubd = they have ^pn = ^ which requires 
a high scalaron mass to comply with solar-system constraints. 

In asymptotically stable scalar-tensor DE models the matter perturbations grow 
slowlier than a in the matter stage and this growth is scale-independent. 

The growth of matter perturbations could provide an efficient way to discriminate 
between modified gravity DE models and DE models inside GR. One can characterize 
the growth of matter perturbations on small redshifts using the quantity 7(2;) 

/ = n^zy^^^ . (58) 

In the pioneering papers on this approach 7 was taken constant [31]. But it is impor- 
tant to realize that 7(2) is a function of z which is generically not constant, and that 
its variation can contain crucial information about the underlying model. As was 
shown earlier, for a wide class of models inside GR one has |7o = ^(-2^ = 0)1 ^0-02 
so that 7(2;) is approximately constant [25j. In ACDM we have 70 ~ 0.55 (with a 
slight dependence on ^^,0) and 7q ^ —0.015. However, 7q can be significantly larger 
in models outside GR [20] . From the growth of matter perturbations we can calculate 
7(2) and ^ and in particular their value at z = 0. As we can see from Figure 2a, we 
find that / is scale-independent on low redshifts z < 0.5. This means that the growth 
of matter perturbations is the same on all relevant scales in this redshift interval. 
This is easily understood as Hq ~ Hq so that all subhorizon scales, and certainly 
the relevant scales that are deep inside the Hubble radius today, satisfy Hq (we 
take ao = 1). Some scale dependence can appear in the growth of matter perturba- 
tions today for scales that are so large that the original equation (H3l) (or (H2l) ) no 
longer provides an accurate approximation. Therefore the quantity 7(2) is essentially 
scale independent too on low redshifts z < 0.3 in our model. On higher redshifts 
some restricted dispersion appears which could be another signature of this model, 
for example aX z = 0.5 we can have a difference A7 ~ 0.04 between various scales, 
see Figure 2b. On even higher redshifts the growth is of course scale-dependent and 
we can see from Figures 2a that the relative increase is damped with increasing n 
in accordance with earlier analytical estimates. As mentioned earlier, this increase 
is constrained by the observations and cannot be large. We find further that in 
this f{R) model 70 = 7(2 = 0) ~ 0.41 which is much lower than in ACDM where 
7o ~ 0.55 (and 7q ~ —0.015). This is an interesting property which clearly allows 
to discriminate this model from ACDM. It is also significantly lower than the value 
found in some scalar-tensor DE models. This value seems essentially independent 
of the model parameter n. So a measurement of 70 could allow to discriminate this 
model from ACDM, but also possibly from other modified gravity DE models. As 70 
is scale independent this is a clear signature of this f{R) model. 

These conclusions remain the same if we let ^2.^,0 vary in the viable range of fim,o 
values. As we can see on Figures 3, there is some variation of 70 and 79 in function 



11 



1.1 

1.0 






0.9 




A=30/7-' Mpc 


0.8 


1 ' 


A=10/r' Mpc 


0.7 


- / / 


A=5/r' Mpc 




/ / 


A=l/r' Mpc 


0.6 


-/ / 

\ 


ACDM 


0.5 


1 

-1 







Figure 2: a) On the left, the growth factor / is shown for the model with n = 2 
and the dimensionless model parameter A has the value A = 0.94. We see the scale 
dependent behaviour of /. On low redshifts all the cosmic scales displayed here 
are in the regime fl5T|) .and matter perturbations obey fl55l) in the matter-dominated 
stage. A significant fifth-force is felt in the deviation of Poisson's equation for matter 
perturbations on these cosmic scales. Note that on higher redshifts these scales will 
eventually be in the opposite regime flS21) . so that perturbations follow the evolution 
( 15^ . This is the origin of the bumps seen on the figure, b) On the right, the quantity 
j{z) is displayed for the same model parameters as on the left panel. We find a low 
value for the growth parameter 70, substantially lower than in ACDM. It is seen that 
the dispersion of 7(2) is very small, while that of 70 is negigible. We find also a value 
for 7q = = 0), 7o ~ 0.2 which is much larger than the value found for DE models 
inside GR. c) On the bottom figure, the quantity is displayed up to 2; = 0.5. 
This quantity does not depend explicitly on the background quantity flm- At 2; = 0.5, 
a difference of about 6% is found with ACDM. 
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0.3/2 


0.315 


0.302 


0.289 


0.273 


0.263 


0.254 


0.245 


0.238 


n 00*7 

0.227 


7o 


0.396 


0.399 


0.404 


0.409 


0.415 


0.419 


0.422 


0.425 


0.428 


0.432 




-0.253 


-0.246 


-0.234 


-0.224 


-0.210 


-0.202 


-0.195 


-0.189 


-0.183 


-0.175 


Za 


0.654 


0.673 


0.711 


0.747 


0.798 


0.830 


0.862 


0.893 


0.922 


0.965 


5^(^=0.5) 

5m, 


0.731 


0.733 


0.737 


0.742 


0.747 


0.751 


0.754 


0.758 


0.761 


0.765 



Table 1: This table summarizes the value of the growth parameters 70, 7o and of z^, 
the redshift when accelerated expansion starts. All these values correspond to n = 2 
and A = 0.94. 



0.44 



0.42 



0.40 



0.38 



0.36 





-0.28 



0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 

n„,.o 



0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 



Figure 3: a) On the left, the parameter 70 is displayed in function of the cosmological 
parameter VLmfl for different values of the model parameter n. For given n, we see 
some variation with the highest value obtained for a low matter density universe. Still 
all these values remain clearly distinguishable from their values in ACDM. b) On the 
right, the parameter 7q. We note that very large values are obtained for all models. 



of Vtmfl- We get 7o ~ 0.4 for Vtmfl = 0.32 and 70 ^ 0.43 for ^^.,0 = 0.23. As for 7q we 
have 7q ~ —0.18 for Qm,o = 0.23 and 7q ~ —0.25 for Qrn,o = 0.32. So a low matter 
density universe brings the quantities 70 and 7q closer to their values in ACDM but 
they still remain far from them. Such low values for 70 seem in trouble with present 
observations (see e.g. fH3]) but an appropriate analysis of the data should relax the 
assumption 7= constant. It is likely that such a low value will remain in tension with 
observations. 

Finally we have here again an example of a modified gravity model where = 
0) = 7g 7^ and actually large. We obtain the high value 7q ~ —0.14 which is largely 
outside the range |7q| < 0.02 found for DE models inside GR. It is also much higher 
than the value found for some scalar-tensor DE models. So we have here again a 
characteristic signature of our model which clearly differentiates it both from ACDM 
and DE models inside GR, but possibly also from other DE models outside GR. Like 
for 7o, the scale independence of 7q makes it a clear signature of the model. We have 
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checked that 7q obeys the constraint 
7^ = [Infi-y 

where we have to include the factor '^"^^^^ ~ ^ 1. As for the scalar-tensor models 
considered earlier, we obtain a nearly linear behaviour on low redshifts z < 0.3. 

In conclusion we find that the growth of matter perturbations on small redshifts 
provides a powerful constraint on our f{R) model. We find low values for the pa- 
rameter 7o with 7o ~ 0.4, and high values for 7q with 7q ^ —0.2. These values are 
not scale dependent and provide therefore a clear signature. The growth parameters 
7o and 7o are mainly affected by the cosmo logical parameter flmfl, in particular 7q is 
less negative for low Qm,o values. But in all cases, characteristic values very far from 
ACDM and from all DE models inside GR are found. 

There are still very large uncertainties on the quantity f{z) or on = ^ where b 
is the bias factor (see e.g. [33], [32]). One should further keep in mind that a precise 
observational determination of both /(z) and flm is needed in order to measure 7(2;) 
accurately. Hence to get precise values for the couple 70, 7o we need to determine 
accurately both f{z) and flm{z) around z = 0. If future surveys will constrain 70, 7o 
to be close to their values for ACDM then the models we have investigated here will be 
ruled out. Though a systematic numerical exploration in the model parameter space 
is very hard to achieve, the model for n = 2 shows already a very large deviation from 
ACDM if one considers the growth of matter perturbations while this model meets 
all other constraints. So either the model will be ruled out for any value of n or the 
growth of matter perturbations will at least significantly restrict the viable interval 
in the model parameter n. 

We conjecture that many if not all viable f{R) models will have similar obser- 
vational signatures. A precise determination of the parameters 70 and 79 could be 
decisive in the quest for the true DE model especially if it is an f{R) modified gravity 
DE model. 
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